Statistical conservation laws in turbulent transport 
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Abstract 

We address the statistical theory of fields that are transported by a turbu- 



Q 
U 

^ ' lent velocity field, both in forced and in unforced (decaying) experiments. We 



propose that with very few provisos on the transporting velocity field, corre- 
lation functions of the transported field in the forced case are dominated by 
statistically preserved structures. In decaying experiments (without forcing 



> 
(N 

o 

Ti^lj- ' the transported fields) we identify infinitely many statistical constants of the 



motion, which are obtained by projecting the decaying correlation functions 
on the statistically preserved functions. We exemplify these ideas and pro- 



^ ' vide numerical evidence using a simple model of turbulent transport. This 

> 

' example is chosen for its lack of Lagrangian structure, to stress the generality 

' of the ideas. 

Turbulent transport represents a dynamical system in which observable quantities (such 
as the concentration of a contaminant) change over time and space in a complex fashion 
|T|]. Unless the transported agent is continuously renewed, irreversible mixing and dispersion 
result in a uniform distribution over the system 0,0]. It had been therefore quite surprising 
that in certain examples (in which the advecting velocity field is (5-correIated in time) it was 
observed that there exists functions of space that remain invariant to the time evolution 
]. In this Letter we propose that with very few provisos on the turbulent velocity field, 
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the dynamics exhibits infinitely many statistically conserved structures. By "statistically 
conserved" we mean structures that are not conserved for every realization of the turbulent 
velocity field, but are conserved after taking an average over all realizations of the turbulent 
velocity field. This observation is very general, pertaining to scalar fields (like contaminants) 
or vector fields (like magnetic fields), with pressure effects or without. Moreover, these 
conserved structures dominate the statistics of the transported fields in the forced case, (in 
which fresh supply is continuously furnished). Indeed, forced experiments are more frequent 
in the turbulence community since they provide much better statistical measurements than 
decaying experiments. The present ideas, if supported by further studies, provide a unifying 
approach to understanding such statistical properties in turbulent transport, including the 
much debated issue of anomalous scaling 0]. 

Consider a turbulent velocity field u{r, t) whose statistics is assumed stationary, but 
without any further restrictions. We consider the wide class of problems in which another 
field is transported passively by that velocity. The field may be an advected scalar 6{r,t), 
with the equation of motion 

^ + u-ve = kvH , (1) 

ot 

or a vector, like a magnetic field B{r,t) satisfying [§] 

dB 



dt 



+ {u-V)B = {B ■V)u + kSJ^B . (2) 



We may also consider advection, as in 0, of a vector field w whose divergence vanishes, 

duo 

— + {u-V)w = -Vp + tiV'^w . (3) 

Other examples can be considered, as long as the equation for the transported field 0(r, t) 
(scalar or vector) has the form 

dt(i) = C(i). (4) 
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Here, £ is a stochastic operator that is built out of the turbulent velocity field. The key 
point is that the equation of motion is linear in the transported field 0. Statistical averages, 
denoted as (...), are performed on the initial conditions and on the realizations of u. A 
fundamental consequence of the linearity of the equations of motion is that the correlation 
functions may be expressed as 

{<f){ri,t)...<f){rN,t)) = 

(5) 

fPS^it) (0(Pi,O)...0(p^,O))orfp, 
where (. . .)o on the RHS is an average over the statistics of the initial conditions. The nota- 
tion r = (ri, . . . , vn) is used for simplicity. Note that we have used the passive nature of the 
transported field, i.e. the fact that the velocity is independent of the initial distribution of 
(j), to separate the averages over the initial conditions and the velocity. The linear operator 
V^^^ (t) propagates the A^th-order correlation function from time zero to time t. 

The evolution operator C generally includes dissipative terms, and without fresh input 
(forcing) the statistics of the field is time-dependent; this is the decaying case (^. A related 
problem of much experimental and theoretical interest is /orcec? turbulent transport where an 
input term / is added to the equation (^. The situations of interest in turbulence typically 
involve an input acting only at large scales of order L. The object of major interest is the 
statistics of the transported field at distances much smaller than L and at the stationary 
state. As usual in turbulent flows, the moments of the velocity increments are assumed to 
behave as power laws with respect to the distance between the two points. The correlation 
functions are then expected to contain anomalous contributions behaving as 

(0(Ari,t) . . .0(Ar;v,t)) = A^^(0(ri,f) . . .0(r^,t)), (6) 

with scaling exponents Civ which cannot be inferred from dimensional analysis. 

The aim of this Letter is to draw attention to two crucial ideas, without further provisos 
on the velocity field. We do not prove these ideas here, we offer them as conjectures, and will 
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exemplify them below using a simple example that strongly differs from the known previous 
cases. Specifically, the model discussed below lacks the Lagrangian structure heavily relied 
upon in previous works. This offers evidence for the generality of the following ideas: 

(i) In the decaying cases, despite the non-stationarity of the statistics, there exist special 
functions Z^^\r) such that 

lW{t)= f Z'^''\r){<p{r^,t)...(^irr,,t))dr (7) 



are statistical integrals of motion. In the limit of infinitely large system it does not change 
with time. In a finite system, and see Fig.l as an example, it is constant in time until some 
outer time scale Tl is reached. It follows from and the conservation of I^^\t) that in 
the infinite size limit the Z^^^'s are statistically conserved structures: 

ZW(I1) = jv^^Xt)Z'^\p)dp. (8) 

(ii) The anomalous part of the stationary correlation functions in the forced problems is 
dominated by statistically conserved structures. A direct consequence is that the small- 
scale statistics of the transported field (j) in the forced case rests on the understanding of 
the decaying problem. A by-product is that the scaling exponents Cjv are universal, i.e. 
independent of the forcing mechanisms for any forcing that is statistically independent of 
the velocity field. 

To appreciate the generality of those statements it is desirable to review briefly the 
development of the ideas leading to them. The first instance in which the role of 

statistically conserved structures emerged was passive scalar advection by a velocity field 
with a short correlation time and self-similar spatial correlations. This is the well known 



Kraichnan model |T^, which is non-generic, but it has the advantage that it lends itself 
to analytic calculations. In the Kraichnan model the simultaneous A^th order correlation 
function satisfies a linear differential equation, which is inhomogeneous in the presence of the 
forcing. The general solution is the sum of the inhomogeneous and the homogeneous parts. It 
turns out that both exhibit scaling properties, but the scaling exponent of the homogeneous 



part (the zero modes) is leading (smaller) compared to the other. Moreover, the subleading 
inhomogeneous exponent is predictable by dimensional analysis a la Kolmogorov, whereas 
the leading exponent is anomalous, and calls for an explicit calculation. In the case of the 
Kraichnan model it is also clear that the homogeneous part of the linear operator governs 
the rate of change of the correlation function in the decaying case, and thus the zero modes 
are statistically conserved structures for the decaying problem. 

Next, the understanding of the nature of the conserved structures followed from the 
analysis of the dynamics of groups of Lagrangian trajectories of tracer particles ||TlHl4|. The 



remark is that for the passive scalar equation ([l|) the transported field is conserved along 
the trajectories of the tracer particles dR{t) = u{R{t),t)dt + \/2ii df3{t) , where /3(t) is a 
Brownian process. To know the scalar field at position r and time t it is enough to track 
the corresponding tracer particle back to initial position p. The evolution operator P^^-* (t) 
in (^ coincides then with the probability density that tracer particles reach the positions 
r at time t given their initial positions p. For example, to understand the exponent Cs 
one needs to focus on the dynamics of three tracer particles. Obviously, three particles 
define at any moment of time a triangle, which in its turn is fully characterized by one 
length scale R (say the sum of the lengths of its sides), two of its internal angles, and all 
the angles that specify the orientation of the triangle in space. When the particles are 
advected by the turbulent velocity field, the scale R of the triangle and its shape (angles) 
change continuously. The statement that can be made is that there exist distributions on the 
space of the triangle configurations, that are statistically invariant to the turbulent dynamics 



rT| , p!2|JT5| . In other words, if we release trios of Lagrangian tracers many times into the 
turbulent fluid, and we choose the distribution of their shapes and sizes correctly, it will 
remain invariant to the turbulent advection. Such statistically conserved structures are the 
aforementioned zero modes and they come to dominate the statistics of the scalar field at 
small scales. The anomalous exponents of the zero modes, such as Ca? can be understood 
as the rescaling exponents characterizing precisely such special distributions. Of course, the 
same ideas apply to any order correlation function with the appropriate shape dynamics. 



The relevance of Lagrangian trajectories can be also demonstrated for the magnetic field 
case (0), by adding a tangent vector to the tracer particle, and see for more details. 

All this was established for the non-generic Kraichnan model of passive scalar advection. 
Generalizing the zero-mode ideas to generic cases has a clear importance, but it calls for 
coping with two different problems. First, if the velocity field in not (5-correlated in time, 
we have no time-independent operators anymore, and the idea of zero modes as functions 
annihilated by a fixed operator cannot be carried over. The correlations functions still obey 
linear equations of motion, but the non-vanishing correlation time of the velocity field sta- 
tistically couples different times. It is therefore not obvious that single-time (simultaneous) 
objects might still be statistically conserved. The first indication in favor of the generality of 
the picture came in [|1^ . The velocity field is generated by the 2-dimensional Navier-Stokes 



equations (in the inverse cascade regime) and has a finite correlation time. It is clear from 
the previous discussion that breaking the non-generic approximation of ^-correlation in time 
is the crucial point. The numerical evidence is that statistically conserved structures are 
still present. Two parallel things have been done: on the one hand, the third order cor- 
relation function (0(ri)0(r2)0(r3)) was measured in forced simulations. This yielded the 
distribution expected to be invariant. On the other hand, trios of Lagrangian trajectories 
were released into the flow, and it was demonstrated that their evolution agreed with the 
notion of statistical invariance. 

The second problem is that the formulations proposed up to now heavily exploit the La- 
grangian structure of the dynamics. In this Letter we take the ideas further and demonstrate 
that their generality transcends the applicability of Lagrangian trajectories. The moment 
that we consider other models of transport, e.g. with pressure, the properties of A^th order 
correlation functions can no longer be connected to N Lagrangian tracer particles. A new 
formulation of the statistical conservation laws is thus needed. That is what (|^ provides 
for. It is clear from it that the key point is not the Lagrangian structure of the correlation 
function propagation but that the linear operator is associated with (statistically) preserved 
functions (cf. Eq.(|^)). The preserved functions then define the integrals of motion. Of 



course, this still does not necessarily mean that they also dominate the field correlations in 
the forced case. To prove it in full generality is difficult, and certainly beyond the scope of 
this Letter. Rather, we will demonstrate the point with a non-trivial example. 

To address the issue in a model as different as possible from those considered so far, we 
analyze shell models, where the discretization of the phase space destroys any Lagrangian 
structure. A possible discretization of the field equation (||) reads |l8l 



dt 



= c 



(9) 



m,m' Ym' 



where the Un variables are generated by the "Sabra" shell model 

du 



dt 



) - uk^Un + fn 



(10) 



where the coefficients a, 6, and c are real. In our simulations k 



5 X 10- 



1, 



b = —0.4, and c = a + b. The wavevectors are A;„ = ko 2™ with n = 0, . . . , N. The smallest 
wavevector is given by ko = 0.05 while defines the ultraviolet cut-off. The velocity forcing 
/„ is limited to the first shell n = 0. For k = u = the energies J2n I'^nP and J2n l^nP 
are dynamically conserved, i.e. realization by realization. The operator p(^) of Eq.(^) takes 
here the explicit form {R{t\to) . . . R{t\to)) where R{t\to) = T+ exp[//^ dsC{s)], with r+ being 
the time ordering operator. 

To demonstrate the statistical conservation laws, two things were done. First we consid- 
ered the forced problem, adding random forcing to the first shell of (P). We have measured 
all the available correlation functions of second, fourth and sixth order. Due to phase sym- 
metry constraints, these are (we put a subscript / to stress that these are statistical averages 
in the stationary forced ensemble): 



F2{n) ^ m')f , 
Fi'\n,m) = (|0„ri0™r)/ 



(11) 
(12) 



Fl >{n) = ((^„+20^+i0^+i'^n-i)/ , (13) 
Fe{n,m,k) = {\(l)n\'\(l)m\'\M')f ■ (14) 

Secondly, we studied the decaying problem, preparing initial states = 0) and following 
their evolution. As initial states we took distributions of 0„ = except for n — 14, 15 where 
we initialized the field with a constant modulus and random phases. We then computed the 
following objects: 

l'-'\t)^j:WmF2{n), (15) 

n 
n,m 

+ E(C+20n+10n+lC_l)(t) Fi'\n) , (16) 
n 

/(6)(t) = Y: {M^I^Mkm) F,{n,m,k) (17) 

n,m,k 

Fig.l summarizes the results. We show, for these three orders, (i) the time dependence of 
the nth order decaying correlation functions themselves, (ii) the time dependence of I^^\t). 
In panel (C) we show also for comparison the time dependence of I^^\t) if we replace the 
measured forced Fq by its dimensional shell dependence (i.e. the shell dependence if the 
Kolmogorov theory were right). We see that only the properly computed /^"^(t) are time 
independent for times smaller than the large scale eddy turn over time, Tl- The decay 
observed for times larger than Tl is simply due to finite size effects intervening when the 
decaying field reaches the largest scales. Indeed, the correlation functions of the forced 
problem (that we used for Z^^^) are not dominated by the inertial contributions at those 
scales. 

In conclusion, the importance of statistical conservation laws for turbulent transport was 
demonstrated. We have exemphfied both the existence of statistically conserved objects and 
the fact that they come to dominate the small-scale behavior of the forced stationary corre- 
lation functions. The only provisos are that the turbulent velocity field has the usual scahng 
properties observed in turbulence, and that the transported field is passive. That offers a 
unifying picture to the statistical theory of turbulent transport. We propose that further 
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research, both numeric and analytic, can contribute to the sohdification of the generahty of 
these ideas. 
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FIG. 1. Panel (A): time dependence of the decaying second order correlation functions (x), 
together with the time dependence of the statistically conserved quantities I*-^-*, (□). Equations 
([lO| ) and (P) have been integrated with a total number of shells N = 33. Time in the horizontal-axis 
is given in units of the eddy turn over time T^. Panel (B): the same as panel (A) but for the fourth 
order correlation function and with = 25. Panel (C): the same of panel (B) but for the sixth 
order correlation function. Here we also present I^^^ when we replace the forced solution FQ{n, m, k) 
with its dimensional prediction, (*) . 
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